Sub-milliKelvin spatial thermometry of a single Doppler cooled ion in a Paul trap 
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We report on observations of thermal motion of a single, Doppler-cooled ion along the axis of 
a linear radio- frequency quadrupole trap. We show that for a harmonic potential the thermal 
occupation of energy levels leads to Gaussian distribution of the ion's axial position. The dependence 
of the spatial thermal spread on the trap potential is used for precise calibration of our imaging 
system's point spread function and sub-milliKelvin thermometry. We employ this technique to 
investigate the laser detuning dependence of the Doppler temperature. 

PACS numbers: 37.10.Vz,37.10.Rs,37.10.Ty 



I. INTRODUCTION 

In the final stages of laser cooling the motion of an 
atom is dominated by the random recoils of photon ab- 
sorption and emission events [jj, l2| . If the atom is har- 
monically confined, this Brownian motion [3|, Ij] is ex- 
pected to result in a Gaussian distribution of its position 
and velocity [sl-Q- The width of this distribution can 
be intuitively interpreted as temperature, which we de- 
fine as the time averaged energy divided by Boltzmann's 
constant for a single particle. This is a powerful no- 
tion, since many experiments require low residual kinetic 
energy, e.g., for precision metrology |8l-[l2| or quantum 
computation and simulation [13|, [ij . In this article we 
study the spatial probability density of a Doppler cooled 
Mg+ ion trapped in a linear radio frequency (rf) trap, 
confirm the expected Gaussian distribution, and demon- 
strate that our straightforward imaging approach enables 
precise thermometry, as required for a wide range of ex- 
periments. 

While in the strong-binding limit the comparison of 
the strengths of motional sidebands allows precise tem- 
perature measurements 115|,[l6|, in the weak-binding limit 
the sidebands are not resolved. In this regime, ion tem- 
peratures are usually derived from the fluorescence line 
shapes which are decomposed into their Lorentzian life- 
time contribution and the thermal distribution [I3 by 
fitting a Voigt function. However, this method relies on 
the assumption of a Gaussian thermal distribution and 
the separation of the lifetime and thermal widths can be 
accompanied by rather large uncertainties |9|, [l8|. Up- 
per limits for the temperatures of cooled ions have also 
been obtained by measuring the (thermal) spatial distri- 
bution in early laser cooling experiments [19] . Uncertain- 
ties down to 5 mK have recently been reported [20] by 
means of a specifically designed Fresnel lens with high 
spatial resolution. A similar technique has also been ap- 
plied to atoms [21|. Thermometry on large ion crystals 
has been performed by comparing crystal images to the 
results of molecular dynamics simulations |22l, |23|. 

In this work (see also [24]), we investigate the time- 
averaged spatial distribution of a single Mg+ ion confined 
in a linear quadrupole trap and laser-cooled close to the 
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FIG. 1: (Color online) Experimental setup: A ^^Mg^ ion is 
trapped in a linear Paul trap. Radial confinement is accom- 
plished by an rf voltage applied to two of the four rods. A 
dc voltage applied to the rings provides tunable axial confine- 
ment. The ion is Doppler cooled by red-detuned laser light 
addressing the cycling D2 transition. A single-photon cam- 
era observes the axial spatial fluorescence distribution via an 
imaging system. 



Doppler limit. The trap is operated with weak axial dc 
confinement which results in an axial spatial spread con- 
siderably larger than the resolution of our imaging optics. 
As expected for Doppler cooling, we observe a Gaussian 
fluorescence distribution. By accurately calibrating both 
the magnification and the resolution of our imaging op- 
tics, we are able to measure accurate values of the ther- 
mal spread which allow to extract ion temperatures with 
sub-mK uncertainties. We employ this precise thermom- 
etry to investigate the laser detuning dependence of the 
Doppler temperature. 



II. THEORETICAL BACKGROUND 

We consider an ion of mass m, trapped along the z- 
axis in a harmonic potential V = muj'^z'^ /2 with secular 
frequency cj, leading to motional quantum mechanical os- 
cillator states tl^n{z) of energy En = (n + l/2)Jiw. Under 
the assumption that the ion's random walk caused by the 
stochastic photon absorption and emission is ergodic, the 
single ion can be assigned a temperature T which quan- 
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FIG. 2: (Color online) Time- averaged spatial distribution 
of a single laser-cooled ion trapped with secular frequency 
uj = 27V X 15 kHz (T ^ 1 mK). The inset shows the ion image 
while the plot is a histogram in axial direction. The residuals 
(lower graph) of a Gaussian fit (solid line) confirm a nor- 
mal distribution as expected from Brownian motion caused 
by the recoils of stochastic photon scattering events. The 
RMS width Azcam of the Gaussian and the RMS resolution 
AzpsF ~ 1 /im of our imaging system are indicated. 



tifies the time- averaged occupation number of the states 
n = {ex.-p[huj / kbT] — 1)~ [25] with the Boltzmann con- 
stant kij. The population probabilities P^ of the states 
V^n(^) follow the distribution P^ = n^ (n + 1)"^''+^^ [25J 
which translates into a time-averaged spatial distribution 
of the ion around its mean position zq 



f{z) = J2Pn\M^ 



(z-zq)^ 



n=0 



27rAzth 



(1) 



In the evaluation of the sum we used Mehler's Hermite 
polynomial formula [26| . The variance of this Gaussian is 
Az^^ = (n + l/2)h/{muj). For k^T :^ hu as appropriate 
for the weak binding regime we have n ~ k^T /(hw) and 
the root mean square (RMS) width 



Azth 



hT 



muj^ 



(2) 



The temperature limit of a Doppler-cooled ion due to 
secular motion in a harmonic potential results from an 
equilibrium of laser-induced cooling and heating rates 
and is given by [23 



^=w^'^^i^'^'^m 



2IAI 



(3) 



It depends on the laser detuning A < and the 
laser intensity / = si sat i where s and I sat are the sat- 
uration parameter and intensity, respectively. F is the 
natural linewidth of the optical dipole transition and 
^ = 2/5 takes the dipole emission pattern into account. 



The temperature diverges for |A| ^ 0, |A| 
s ^ oo. The minimum of Tmin = \/l + 5(1 
is obtained at a detuning of Amin = — FyT^ 
for small 5, reduces to Amin ^ —T/2. 



-^ oo, and 
±^)hr/4kj, 
- s/2 which, 



III. EXPERIMENTAL SETUP 

Our experimental setup is shown in Fig. [1] (see 
also [11|). A single ^^Mg+ ion is trapped in a linear 
rf quadrupole trap which operates at a trap frequency 
1] = 27r X 22.6 MHz and generates radial rf confinement 
with a secular frequency cj^ ^ 27r x 1 MHz. The trap 
electrodes are surrounded by two rings with a dc volt- 
age applied to generate axial confinement tunable from 
cj ~ 27r X 10 to 150 kHz. These frequencies are measured 
by secular excitation with a weak external signal and can 
thus be controlled with an accuracy of Au < 27r x 1 kHz. 
Since Ur '> uj^ axial and radial motion are decoupled 
so that we can neglect radial movement in the follow- 
ing. Due to its zero nuclear spin, the alkali-like spec- 
trum of ^^Mg+ shows no hyperfine structure and the 
Di and D2 lines constitute clean cycling transitions well 
suited for Doppler cooling without the need for repumper 
lasers. Two laser beams each stabilized in frequency and 
intensity address the 3^5'i/2-3^P3/2 D2 transition near 
280 nm (natural linewidth F = 27r x 41.8(4) MHz 0, 
I sat = 2.50 kW/m^). One beam is aligned along the 
axial trap direction. The second laser beam is slightly 
tilted by 14° against the first beam to provide radial 
cooling. The beams are detuned with respect to each 
other by ~ 500 kHz to avoid a stable interference pat- 
tern. The ion is imaged with a //2 condenser lens and a 
microscope objective onto a single-photon camera (SPG, 
Quantar Mepsicron II). Because of the SPC's limited spa- 
tial resolution of 56 /im, a magnification of M ^ 100 is 
chosen for the imaging system. The detector plane is 
digitized in pixels of 49 /im size. The conversion fac- 
tor between the real space object size in /im and image 
size in pixels is calibrated accurately by measuring the 
distance between two simultaneously trapped ions in the 
imaging plane of the camera for several trapping poten- 
tials uj [29]. The RMS resolution of the imaging system's 
point spread function (PSF) at 280 nm is Azpsf ^ 1 /^m 
for optimal alignment. It derives from the resolution of 
the lens system, the resolution of the camera, and the 
discretization of the camera data. However it is dimin- 
ished when the ion is shifted out of focus, particularly by 
varying laser forces, and thus becomes slightly detuning 
dependent. After conversion into /im in object space, the 
RMS spot size Azcam of a recorded ion image appears as 
a convolution of the thermal spread Azth with the finite 
PSF width AzpsF' Assuming the PSF to be Gaussian, 
we have: 
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FIG. 3: (Color online) (a) Measurement of the spot size Azcam of the ion images with A = — 27r x 18.7 MHz and s ^ 0.1 for 
various trapping potentials uj /2ti. The insets show exemplary ion images. From a fit (solid line) of Eqs.[2]and[4l we obtain an 
ion temperature of T = 1.02(3) mK and a PSF of Azpsf = 1.13(3) /xm RMS with statistical uncertainties, (b) Same analysis 
for various laser detunings A. Results are shown as squares together with temperatures expected from Eq. [3] (dashed gray line). 
We attribute the observed temperature excess of up to ^ 0.2 mK to systematic micromotion, see discussion in the text. The 
values of AzpsF obtained for all data points agree within their statistical uncertainties. This measurement demonstrates the 
precision of the thermal spread thermometry providing a total accuracy of < 0.3 mK unchallenged by other methods in the 
unresolved sideband regime. 



IV. MEASUREMENT AND RESULTS 
A. Time-averaged spatial distribution 

The inset of Fig. [2] shows an image of a Mg+ ion 
trapped with a secular frequency a; = 27r x 15 kHz and 
cooled close to the Doppler limit of about T ^ 1 mK 
(n ^ 1400). From Eq. E] follows a RMS spatial spread 
of Azth ^ 8 A^m which is about eight times larger than 
the PSF of our imaging system. In this case, with 
Eq. m the PSF only contributes ^ 1% to the width 
^Zcam of the image. The laser detuning is set to 
A = -27r X 40 MHz ^ -F, and the total laser intensity 
- with both laser intensities equal at the location of the 
ion - is limited to s ^ 0.1 which is monitored by the 
observed photon scattering rate. Figure [2] shows a his- 
togram of the axial spatial distribution obtained from 
the image. No statistically significant deviation from a 
Gaussian could be found with a fit, as can be seen from 
the residuals, which confirms the stochastic nature of the 
ion's motion. The RMS spot size of the ion image of 
^Zcam = 7.7(1) /im is indicated by the horizontal bar, 
while the vertical bar shows the resolution of our imag- 
ing system. 



(Eq. [2]). From measurements of the spot sizes Azcam for 
different cj, we obtain Azpsf-, ^Zth^ and thus the abso- 
lute temperature T of the ion from a fit of Eqs. [2] andJH 
This method depends on the shape of the PSF and the 
constancy of its width Azpsf- The ion images at high 
secular frequencies which reflect the PSF do not show a 
significant deviation from a Gaussian. To ensure the con- 
stancy of its width, we readjust the imaging system for 
optimal resolution at each data point, thus compensating 
the effects of ion position changes due to the varying bal- 
ance between the trap potentials and the laser force. At 
the same time we minimize radial micromotion at each 
data point to avoid axial heating via possible coupling to 
the radial motion. Figure [3] (a) shows such a measure- 



ment for A 



-27T X 18.7 MHz and s < 0.1 with uj be- 



tween 27r X 12 and 124 kHz. The insets show correspond- 
ing images and the resulting axial RMS spot sizes Azcam 
obtained from Gaussian fits to the image histograms are 
plotted versus uj/27r. Note that the radial widths only in- 
crease slightly by about 15 % towards lower axial poten- 
tials because radial micromotion compensation becomes 
increasingly difficult. The 1/uj behavior of Azth is con- 
firmed by a fit (Eq. [3j), which yields an ion temperature 
of T = 1.02(3) mK and a PSF of Azpsf = 1.13(3) /im, 
both with statistical uncertainties. 



B. Spatial thermometry measurement 

In order to obtain a precise value of the absolute tem- 
perature for larger oscillation frequencies uj as well, we 
need to determine Azpsf niore accurately. We use the 
fact that, according to the laws of Brownian motion, 
the width Azth of the thermal distribution varies ex l/o; 



C. Systematic effects 

The main systematic uncertainty arises from resid- 
ual axial micromotion caused by the axial component 
of the rf fringe fields. The simultaneous presence of 
dc and rf fields leads to a parametric potential V{t) = 



^mft'^{a — 2q cos{flt)) z'^ . ft is the trap frequency, a ex 
elide /'^^'^ and q oc eVrf/mQ'^ are trap parameters as- 
sociated with the dc and rf voltages Udc and K-/, re- 
spectively. Without damping and diffusion, the ion's 
equation of motion is a Mathieu equation and for cer- 
tain parameter sets (a, q) there are stable trajectories of 
the ion z{t) = ZsCos{ujt){l + ^cosQt) with amplitude 
Zs and uj ^ ^Q{a + q^/2y^'^ [30]. The latter can be in- 
terpreted as the secular frequency of the ion in a time 
averaged pseudopotential. For vanishing a it reduces to 
pure rf confinement ujrf = q^/VS while pure dc confine- 
ment uodc = >/a^/2 results for g' = 0. Superimposed on 
the secular motion, the ion performs micromotion at the 
trap frequency Q with an amplitude z^ = qZs/2. For an 
ion cooled close to the Doppler limit of 1 mK, z^ is of 
the order of 50 nm and does not influence the apparent 
spread or lead to significant deviations from a Gaussian 
distribution |6l. l31j. 

A stronger effect due to micromotion is expected via 
the kinetic energy which is given by ^30i] 



Eki: 



,2^2 



-mu z^ + ^mcj^/^5 



zf + 



y-^Azl (5) 



The first term is due to secular motion, the second re- 
flects the contribution from the unavoidable ordinary mi- 
cromotion that is associated with the secular motion. It 
depends on the secular amplitude Zs and is significant 
for low secular frequencies u ^ uJrf- The third term is 
caused by a mismatch Azq between the rf and dc po- 
tential minima positions. It results in a displacement 
5z = AzQCUdc/'^ of the ion out of the rf minimum so that it 
is exposed to excess micromotion [30| that contributes an 
energy E^^^^^^ = ^muj'^rdz'^. A cooling laser introduces 
both damping and a stochastic force to the equation of 
motion. Blatt et al. ^] have solved the corresponding 
Fokker-Planck equation for the first two terms in eq. [5l 
which results in Gaussian spatial and velocity distribu- 
tions. It turns out that the presence of a stochastic force 
transforms part of the micromotion energy into irregular 
motion (temperature) by rf heating which is reflected in a 
broadening of the velocity distribution as well as the spa- 
tial distribution as compared to the case of pure secular 
motion. 

From the lowest secular frequency that we have ob- 
served when lowering the dc potential a, we estimate 
an upper limit for the residual axial rf confinement of 
UJrf = 27r X 7 kHz, corresponding to q = 8.8 x 10~^. 
According to [6], the corresponding ordinary rf heating 
broadens the spatial distribution by a factor of 1.15 for 
the lowest secular frequency of cj = 27r x 12 kHz. This 
broadening decreases strongly towards higher uj :^ ujrf- 
For the estimated upper limit of uJrf = 27r x 7 kHz, our 
Eq. [2] thus overestimates the temperature by a factor of 
1.3 which results in a systematic uncertainty of ^ 0.3 mK 
in our experiment. 
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FIG. 4: Measured ion image spot size Azcam versus laser de- 
tuning for s ^ 0.1 and fixed cj = 27r x 32 kHz. The black 
solid line shows the result expected from Doppler cooling ac- 
cording to Eqs. O together with Eq. [2] and |4] using s = 0.1 
and AzpsF = 1-5 /xm. Qualitative agreement with the over- 
all temperature behavior of Eq. [3] is found. However quan- 
titatively the measured spreads and thus the temperatures 
are about 10% higher than expected which we attribute to 
rf heating due to excess micromotion, which depends on the 
photon scattering rate and thus decreases towards larger de- 
tunings. The gray long-dashed line represents Doppler cool- 
ing with an additional detuning-independent heating rate, the 
short-dashed line assumes a broadening of the line to 1.1 F 
as caused, e.g., by magnetic fields. Both assumptions can not 
explain the observed thermal spread. Larger error bars for 
larger negative detunings are due to lower photon count rate. 



D. Thermometry of a Doppler cooled ion 

In the next step, we employ our thermometry method 
to investigate the detuning dependence of Doppler cool- 
ing. Figure [3] (b) shows absolute temperatures measured 
for laser detunings A between —0.2 and — 1.0 F indicated 
by squares while the circles are the Doppler temperatures 
predicted for pure secular motion by Eq. [3l Our measure- 
ments follow the detuning dependence of Eq. [3] qualita- 
tively, but show an offset of ~ 0.2 mK. The higher tem- 
peratures are on the order of our conservative estimate for 
the systematic uncertainty. However, such a temperature 
rise can also be caused by excess micromotion in a signif- 
icantly smaller residual axial rf potential [30| . While the 
ordinary micromotion becomes insignificant for cu ^ cj^/, 
excess micromotion due to a large A^^o ^ Zg leads to a uj- 



independent contribution E^f^^^^ ^ ^mou'^rAzQ in Eq. [5] 
and may thus increase the ion temperature through rf 
heating correspondingly. 

In order to study the temperature over a wider range 
of detunings between and —3.5 F, we measure the 
image spot sizes Azcam of the ions for fixed values of 
cj = 27r X 32 kHz and s ^ 0.1 (Fig. g]). At this secular 
frequency the ordinary micromotion increases the ion 
temperature by less than 1%. Due to possible variations 
of the resolution Azpsf during the laser scan, we 
do not extract the temperature in this measurement. 



Nevertheless, the spot size A^cam can stih serve as a 
decent measure for the temperature. The sohd hue 
in Fig. m represents the behavior according to Eq. [3] 
in combination with Eqs. [2] and HI using 5 = 0.1 and 
^zpsF = 1.5 jivd. The latter is estimated to be an upper 
limit. Quantitatively, the measured values are again up 
to 10% higher. The data supports our hypotheses that 
rf heating due to excess micromotion is the main cause 
of the higher observed temperatures. Since rf heatingis 
related to the photon scattering of the cooling laser [6| , 
its rate is expected to be detuning-dependent similar 
to the cooling rate. The wide tuning range in this 
measurement allows to exclude two other possible effects 
as causes for the excess temperature. As a comparison, 
the long-dashed line shows the expected broadening 
taking an additional, detuning-independent heating 
rate into account. The short-dashed line assumes that 
the line width is homogeneously broadened to 1.1 T 
caused, e. g., by the Zeeman splitting or micromotion- 
induced line broadening [30|. From earlier spectroscopy 
experiments [9|, we know that for our experimental 
parameters the homogeneous broadening is in fact lower. 
In conclusion, both effects can be excluded as causes for 



the observed higher thermal spreads (temperatures). 



V. SUMMARY 

To summarize, we have demonstrated that a single, 
Doppler cooled ion weakly bound in an ion trap, shows 
a Gaussian thermal spatial distribution as expected from 
Brownian motion. We show that this time-averaged spa- 
tial distribution can be used for sub-milliKelvin ther- 
mometry with an accuracy of < 0.3 mK unchallenged by 
similar methods in the unresolved sideband regime. Fur- 
ther, we have employed this thermometry to investigate 
the detuning dependence of the temperature in Doppler 
cooling. Note that this method which was demonstrated 
for the axial trapping direction here, can in principle be 
applied to any projection of an ion on the imaging plane. 
The use of high NA objectives with resolutions close to 
the diffraction limit [20|, l32| would allow to further im- 
prove on accuracy and to perform thermometry in steeper 
trapping potentials. 
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